The interaction of gravitational waves with strongly magnetized plasmas 
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We study the interaction of a gravitational wave (GW) with a plasma that is strongly magnetized. 
The GW is considered a small disturbance, and the plasma is modeled by the general relativistic 
analogue of the induction equation of ideal MHD and the single fluid equations. The equations are 
derived without neglecting any of the non-linear interaction terms, and the non-linear equations are 
integrated numerically. We find that for strong magnetic fields of the order of f 15 G the GW excites 
electromagnetic plasma waves very close to the magnetosonic mode. The magnetic and electric field 
oscillations have very high amplitude, and a large amount of energy is absorbed from the GW by the 
electromagnetic oscillations, of the order of fO 23 erg/cm 3 in the case presented here. The absorbed 
energy is proportional to Bq, with Bo the background magnetic field. The energization of the 
plasma takes place on fast time scales of the order of milliseconds. The amount of absorbed energy 
is comparable to the energies emitted in the most energetic astrophysical events, such as giant flares 
on magnetars and possibly even short Gamma ray bursts (GRB), for which the mechanism analyzed 
here also has the fast time-scales required. 
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Gravitational waves (GW) can carry a large amount 
of energy near the sources where they are generated (e.g. 
|10|). They tend not to interact much with matter un- 
der normal conditions, it has been shown though in a 
number of articles (e.g. ||, Q, 0, 0, 0, [3, 0, 
[3 0) that GWs excite various kinds of plasma waves, 
the more efficient, the stronger the background magnetic 
field is and the more tenuous the plasma is. Most of 
these studies are analytical and the equations describing 
the GW-plasma interaction were linearized, only Ref. M 
made an analytical study of a non-linear model, Refs. |lj 
and j9l took some second order effects into account, and 
Ref. |3 performed a numerical study. The GW-plasma 
interaction is a totally non- linear effect, and there is so- 
far no conclusive answer to the question of how much 
energy can be absorbed by a plasma from a GW. 

Here, we study the GW-plasma interaction in its full 
non-linearity, solving the non-linear system of equations 
numerically. Our main interest is in the amount of energy 
absorbed by the plasma from the GW and in the kind of 
plasma waves excited by the GW, and we focus on the 
case of very strong magnetic fields of the order of 10 15 G. 

The interaction of gravitational waves (GWs) with 
plasmas, if it is efficient in transferring energy from the 
GW to the plasma, might well be the basic mechanism 
behind the most energetic astrophysical events, such as 
giant flares on magnetars (i.e. highly energetic outbursts 
in the stars' magnetospheres; see e.g. [2(|) or possibly 
even short Gamma ray bursts (GRB; see e.g. [T||; mag- 
netars are strongly magnetized neutron stars [see e.g. jj|], 
which at the same time appear as Soft Gamma Ray re- 
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peaters [SGR]). The origin of short GRBs is far less es- 
tablished than that of the long GRBs. One model re- 
cently discussed is that short GRB might be giant flares 
on magnetars. In this model, the flares are caused by a 
catastrophic reconfiguration of the stellar magnetic field 
(e.g. 0, 0, 0). The energy released in short GRBs 
is currently estimated to be at least f0 48 erg 0], the 
uncertainty is mainly due to the problem of associating 
short GRBs with physical objects. The duration of short 
GRBs is less than 2 seconds. Giant flares have similar du- 
rations and release energies in the range 10 44 to f O 46 erg 
(e.g. 0). 

Our results show that a GW generated by a magnetar 
might well be the source of the energy released in a gi- 
ant flare and may-be even in a short GRB, through the 
absorption of the GW energy by a plasma in the vicin- 
ity of the magnetar. In this scenario, the GW-plasma 
interaction is the primary mechanism that energizes the 
plasma before a giant flare or may-be a short GRB. Giant 
flares and possibly short GRBs would then be an indirect 
observation of GWs. 

The GW is considered as a small amplitude perturba- 
tion of the otherwise flat spacetime, and we assume it to 
be + polarized and to propagate along the z-direction, so 
that the metric has the form g a b = diag(— f , 1+h, l—h, 1), 
with h(z, t) << 1 the amplitude of the GW [6j. Our aim 
is to express the final equations in terms of the potentially 
observable quantities (electric field E, magnetic field B, 
3- velocity V of the fluid, and rest-mass density p) , which 
can either be defined in a Local Incrtial Frame (LIF) or in 
an orthonormal frame (ONF) |6j. Here, we use the ONF, 
because it is a global frame, and it can be shown that the 
ONF in our case is locally equivalent to a LIF when ap- 
plying the particular coordinate transformation given in 
[ll| . The ONF has the advantage that 4- vectors and ten- 



2 



sors take the same form as in flat space-time, the effects 
of curvature appear only through the covariant deriva- 
tives, which are calculated with use of the Ricci rotation 
coefficients [TlJ . Indices of quantities in the ONF carry 
a hat in the following. 

We assume an ideal conducting fluid, so that the 
electric field is given by the ideal Ohm's law, = 
F ah u~ h jc, which in the ONF takes the usual form, = 



7LE + \Vx Bj, with 7 = 1/yi - V 2 /c 2 , F &b Fara- 
day's field tensor, u a the 4- velocity, u a = 7(0, V x ,V y ,V z ), 
and c the speed of light. The evolution of the mag- 
netic field is determined b y t he Maxwell's equation, 
F ah£ + F bc-A + F ca;h = « Gil • electromagnetic 
(EM) energy momentum tensor is defined as T^ M ^ = 

F b c ~ h^F^Fcd) > and for the nuid > we have 
the energy momentum tensor = Hu a u b + rj ah pc 2 , 

where H is the enthalpy, p the pressure, and = 
diag(— 1, 1, 1, 1) the metric of flat spacetime We as- 
sume an ideal and adiabatic fluid, so that H — pc 2 + 
+ p, with r the adiabatic index. The total energy 
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momentum tensor T 
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yields the mo- 
mentum and energy equations T ab - = Continuity 
is expressed by (pu a ) .. = 0. The evolution of the GW is 
determined by the linearized Einstein equation, where we 
take the back-reaction of the plasma onto the GW into 
account 
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where ST XX , ST yy are the non-background, fluctuating 
parts of the components T xx , T yy of the total energy mo- 
mentum tensor, and G is the gravitational constant 0. 
To close the system of equations, we assume an adia- 
batic and isentropic equation of state, p = K p T , with K 
a constant. 

We focus on the excitation of MHD modes which prop- 
agate in the z-direction, parallel to the propagation di- 
rection of the GW and perpendicular to the background 
magnetic field Bq = B$e x . We let consequently B\\e x 
E\\e y and V^e^, and all variables depend spatially only 
on z (note that B x in the following is the total magnetic 
field, it includes Bq). In specifying the general equa- 
tions to this particular geometry, (i) we express all 4- 
vector and tensor components through the potentially 
observable B x , V z , and p; (ii) we expand the covariant 
derivatives; (iii) we keep all non-linear terms, no approx- 
imations are thus made (except for the linearized Ein- 
stein equation). In this way, we are led to a system 
of non-linear, coupled, partial differential equations in 
a spatially 1-D geometry: With the electric field E y from 
Ohm's law, E y = —^V Z B X , Faraday's equation is fully 
expanded to 



d t B x = cd z E„ 
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Expansion of the ^-component of the momentum equa- 
tion yields 



d t q z + d z 
+d z 



4tt 



{-E y B x ))V 2 



8^ + ^ 



- c 2 d zP 

d h 

~B X (cB x + E y V z ) JY^h) 
~E y{c E y + B x V z)j ^- 



-q z (d t h + V z d z h) 



h 



l-h 2 



0, 



(3) 



where we defined the new momentum variable q z as q z := 
HVzj 2 + (—E y B x ). The continuity equation takes the 
form 



dtD + d z (DV z ) - (Dd t h + DV z d z h) 



l-h 2 



0, (4) 



with the new density variable D := -yp. The GW evolves 
according to 

dak = c 2 d zz h+ ^ (E 2 {Bl B 2 )) + ^(p- P0 )h, 

where the background magnetic field Bq and background 
pressure po = Kp^ have been subtracted. 

To recover 7, p, p, H, and V z from the explicitly evolv- 
ing variables q z , D, and B yi we solve the definition of q z 
for V z and insert it into a reformulated definition of 7, 
which yields a non-linear equation for 7, which we solve 
numerically. Once 7 is recovered, all the other primary 
variables follow in a straightforward way. 

We solve the GW-plasma system of equations apply- 
ing a pseudo-spectral method that is based on Cheby- 
shev polynomials (see e.g. 0)- Time stepping is done 
with the method of lines, using a fourth order Runge- 
Kutta method with adaptive step-size control. The one- 
dimensional grid along the z-direction consists of 256 
grid-points and corresponds to a physical domain along 
the z-axis of length L = 5.4I0 7 cm. The sampling time 
step At is set to At = T gw /1A, with T gw — 1/ f gw and 
f gw the GW frequency. 

We assume a background magnetic field B of I0 15 G, a 
background density po = 1CP 14 g cm~ 3 , and an adiabatic 
index T = 4/3. The initial conditions are B x (z,0) = Bq, 
V z (z,0) = 0, p(z,0) = po, and h(z,0) = 0. The GW 
has as boundary condition at the left end zl of the box 
h{zL,t) = /io(i) cos(k gw ZL — ujgivt), so that a monochro- 
matic plane wave is entering the box. The amplitude 
Hq rises within roughly 1 ms from to KT 4 , at which 
value it stays constant for about 3 ms, where after it 
decays to again. At the right edge z« of the box, 
we apply non-reflecting boundary conditions to h(zR,t). 
B x , V z , and p have free outflow boundary conditions at 
both edges of the box. The GW frequency is set to 
fgw — 5 kHz, and the GW dispersion relation is of the 
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FIG. 1: Top: Average energy density e gw (t) of the GW as 
a function of time. — Bottom: Mean total energy density 
e P i (t) of the plasma as a function of time. 



form 2nfg W = u gw — k gw c, with k gw the wave- number of 
the GW. 

The total mean energy density e p i in the plasma at a 
given time t is numerically determined as 



e p i(t) = 



±- j E v { Zl t) 2 dz + ±- I {B x (z,t) - B ) 2 dz 
/L (6) 



+| / p(z,t)v z (z,t) 2 dz 



with L the size of the system — note that we sub- 
tract the constant background magnetic field Bq in or- 
der to take into account only the energy that is in the 
wave motion. Fig. ^ shows e p i (t) and the mean energy 
density e gw (t) of the GW as a function of time, where 

e gw{t) = -^Q0j gw h(t) 2 , with h(t) the mean instanta- 
neous amplitude of the GW oscillation, defined as the 
root mean square average over the entire simulation box 
0. At maximum GW amplitude, the energy density of 
the GW amounts to e gw — 1.3 x 10 27 erg/cm 3 . Once 
the GW enters the system, the plasma starts to absorb 
energy from the GW, and in roughly 2 ms, slightly de- 
layed in the beginning but finally in parallel with the GW 
reaching its maximum energy, the absorption has reached 
its maximum, the energy density in the plasma is roughly 
10 23 erg/cm 3 . The absorbed energy is a fraction 10~ 4 of 
the GW energy density, so that the back-reaction onto 
the GW is not yet important. When the GW leaves the 
system, the energy in the plasma decays almost together 
with the GW amplitude, i.e. the excited waves propagate 
out of the simulation box. 



FIG. 2: Spatial Fourier transform \E v (k z , t)\ 2 of E y (z,t) at 
a fixed time t = 0.00545 s. 



The GW excites wave motions in the plasma that 
travel with the GW and whose amplitudes increase lin- 
early towards the out-flowing edge of the box. Maxi- 
mum amplitudes attained are 3 x 10 12 statvolt/cm (9 x 
10 16 V/m) for the electric field, 3 x 10 12 G for the oscillat- 
ing magnetic field, roughly three orders of magnitude less 
than Bq, and finally the fluid velocity oscillations have a 
maximum amplitude 8 x 10 7 cm/sec. The fluid motions 
thus remain non-relativistic, so that our non-relativistic 
estimate of the kinetic energy is justified. 

The waves excited in E y , B x , and V z have wave- 
number k z and frequency lu that cannot be distinguished 
from k gw and u) gw within the numerical precision of the 
simulation, see Fig. [21 In particular, we do not find any 
harmonics to be excited. The relativistic Alfven speed 
u A = B 2 / '(4-7T po + Bq/c 2 ) is very close to the speed of 
light ((c — ua)/c = 2 x 10~ 24 ) so that the excited plasma 
modes are indistinguishable from magneto-sonic modes. 

The numerical results show that, for strong magnetic 
fields, a large amount of energy is absorbed by the plasma 
from the GW on a short time-scale, which is of the order 
of 10 GW periods, i.e. in the millisecond range. The en- 
ergy absorbed by the plasma is thus not proportional to 
the duration of the GW-plasma interaction if this dura- 
tion is longer than typically a few GW periods. 

In a parametric study, we found that the energy ab- 
sorbed by the plasma is proportional to B 2 (see Fig. |3J) 
and to h 2 ,. Varying p in the range 10~ 20 g/cm 3 < p < 
10 5 g/cm 3 , it turned out that the absorbed energy is inde- 
pendent of the value of the matter density po, see Fig.0J 
the kinetic energy is actually negligible compared to the 
electromagnetic energy (this is in accordance with the 
fact that the background matter rest-energy density is 
much smaller than the background magnetic energy den- 
sity in the entire range of values po investigated). The 
absorbed energy density is furthermore proportional to 



as we verified by varying the frequency in the kHz 
range (from 1 to 10 kHz). It also seems that the time for 
the plasma needed to reach the maximum level of energy 
absorption is related to the time needed for the GW to 
cross the box, L/c, which equals 0.002 sec for the case 
considered here. 

The maximum amplitudes of the excited oscillations 
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FIG. 3: The energy density e p i absorbed by the plasma (at 
maximum absorption) as a function of the background mag- 
netic field i?o (diamonds), together with a reference line of 
logarithmic slope 2 (dotted). 
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FIG. 4: The energy density e p i absorbed by the plasma (at 
maximum absorption) as a function of the background matter 
density po (diamonds, connected with a dotted line). 

are proportional to the box length L, as explained above. 
This implies that the amount of energy density absorbed 
is proportional to the squared box- length, L 2 , with the 
physical meaning of L to be the length along the propa- 
gation direction of the GW where the GW meets a con- 
stant magnetic field. We additionally verified this scaling 
behavior with numerical simulations in which the box- 
length L was varied. The absorbed energy is thus propor- 
tional to L 3 and to the effective area A e f f through which 
the GW is incident on the plasma, where A e ff := V/L, 
with V the volume in which the interaction takes place. 
We can summarize our numerical findings for the to- 
tal energy E p i absorbed by the plasma as follows, not- 
ing that in the case presented here the absorbed energy 
through an effective area of 1 cm 2 is L x 10 23 erg, with 
L = 5.4 x 10 7 cm, 

e ^ 3 -- 1I)7 (tI;) 3 (t^)(toSg) 2 

and E p i is independent of po ■ 

Astrophysical Application. On the basis of our results, 
we can suggest a new model for the primary mecha- 
nism behind giant flares on magnetars and possibly even 



FIG. 5: Sketch: The GW interacts with a plasma volume in 
the equatorial plane (see text for details). 



1 c 50 




1 10 
r,Jr. 



FIG. 6: Estimate of the energy E p i absorbed by the plasma 
as a function of the inner distance r» n from the star, for dif- 
ferent interaction lengths L (long dashes: L = 100 km, short 
dashes: L = 250 km, solid: L — 500 km, and dash-dotted: 
L = 1000km), and for a = 10°, (p = 180°. 



behind short GRBs. This mechanism deposits a large 
amount of energy in a plasma on time-scales of milli- 
seconds, prior to the actual outburst of the giant flare or 
possibly the short GRB. In our model, we assume that 
a magnetar generates a GW, which travels away from 
the magnetar and interacts with a plasma in the vicinity 
of the magnetar, whereby the plasma efficiently absorbs 
energy from the GW. Once the plasma is energized, sec- 
ondary mechanisms will be triggered that convert the 
energy ultimately to Gamma rays, the nature of these 
mechanisms is though not in the scope of this article. 

The typical magnetic field of a magnetar is cur- 
rently estimated to be of the order of 10 15 to 10 16 G at 
the surface (e.g. 0, 0, HO])- GWs are expected to be 
emitted by magnetars, as by usual neutron stars, when 
they undergo some deformation, e.g. as the result of the 
rearrangement of the strong stellar magnetic field. In 
such a case, the GW amplitude at 5 stellar radii typi- 
cally is 10 -4 (e.g. 0| )• At the same distance, the back- 
ground spacetime can be considered flat. The magnetic 
field is poloidal within the light cylinder, so that, near the 
equatorial plane, a GW generated by the magnetar and 
traveling radially outwards propagates perpendicular to 
the magnetic field, as in the presented simulations. 

In order to get a rough estimate of the energy absorbed 
by a plasma, we consider the plasma in the conical vol- 
ume of radial size L between the inner and outer radii r^ n 
and r ou t = ri n +L, respectively (with r the distance from 
the magnetar), which is limited to within the poloidal 
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opening angle a and the toroidal range cp, as illustrated 
by the sketch in Fig. [5] We assume the magnetic field 
and the GW amplitude to be constant over the length 
L, with values B = B^r* /r in ) 3 and ho = Kfa / r in ) , 
respectively, in order to be able to apply Eq. J7J (note 
that the magnitude of the magnetic field along the radial 
direction is proportional to 1/r 3 ). The plasma volume is 
given as V(r m ,L) = § <psm(a/2) (r 3 ut - rf n ), with cor- 
responding effective area A e ff = V(r in ,L)/L. 

For a numerical estimate, we set B+ = 10 16 G, h* = 
1CP 4 , r* = 10km, a = 10°, and, to take possible 
anisotropic GW emission into account, ip = 180°. Fig. 
shows the absorbed energy as a function of rj„ for four 
different values of L. The energy decreases fast with 
distance, it roughly is proportional to l/rf„, so that at 
Tin — 5 r *j the energy has fallen to E p i = 10 43 erg for 
e.g. L = 500 km. Extrapolating our results to regions 
closer to the star, we find at ri n = 2r+ and again for 
L = 500 km a plasma energy E p i — 10 45 erg, which is of 
the order of the energy released in giant flares on mag- 
netars (e.g. Q, 0, |2(j). Even closer to the star, at 
Tin = r *> the energy is of the order of E p i — 10 47 erg, 
well in excess of the energy in giant flares, and approach- 
ing the energy observed in short GRBs according to latest 
estimates (e.g. ^(J, 0]). For different values of B+ and 
h*, the energy values given here scale according to Eq. 
l(7)l in a straightforward way. 

Again for L = 500 km, the volume is 2 x 10 7 km 3 
for ri n = r* and increases to 4 x 10 7 km 3 for ri n = 
lOr*, which corresponds to effective areas of A e tf = 
5 x 10 4 km 2 and A e ff = 8 x 10 4 km 2 , respectively. The 
involved volumes and areas are thus relatively small. 

Our numerical simulations were done for the case of 
a flat background spacetime, which does not hold in the 
range r% n < 2r_, t anymore, to which we have extrapolated 
the energy estimates. The non-flatness of spacetime close 
to the star, together with the necessarily arbitrary choice 



of a value for L, imply that our energy estimates must 
be interpreted with some care, they can though be taken 
indicative of the fact that the GW-plasma interaction is 
efficient and is an important mechanism near the star. 
Also to mention is the uncertainty concerning the actual 
magnetic topology in the flaring magnetosphere, since 
any deviation from the dipole is likely to intensify the 
GW-plasma interaction through the enlargening of the 
possible interaction regions. It thus remains to be seen 
in how far the given numbers will be modified when a 
more realistic decay of the background magnetic field is 
used in the simulations and when the curvature of the 
background spacetime is included, which though only can 
be done when extending the equations to the case of three 
spatial dimensions. In favor of the model of giant flares 
on magnetars (and, less certain, but possibly also of short 
GRBs) driven by GWs is that the mechanism proposed 
has a fast enough time-scale, of the order of milli-seconds. 

In summary, our results show that strongly magne- 
tized plasmas are efficient absorbers of GWs, largely ir- 
respective of the plasma density, and with an absorption 
time-scale of the order of milli-seconds. This implies that 
GWs may be the energy source for secondary, highly en- 
ergetic phenomena. It also implies that GWs may even- 
tually be strongly damped, if appropriate conditions are 
met. In particular, we can conclude that the GW-plasma 
interaction is an efficient and important mechanism in 
magnetar atmospheres close to the star. 
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